Crystal-fields in YbInNi 4 determined with magnetic form factor 
and inelastic neutron scattering 
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The magnetic form factor of YblnNu has been determined via the flipping ratios R with polarized 
neutron diffraction and the scattering function S(Q,w) was measured in an inelastic neutron scat- 
tering experiment. Both experiments were performed with the aim to determine the crystal-field 
scheme. The magnetic form factor clearly excludes the possibility of a TV doublet as the ground 
state. The inelastic neutron data exhibit two, almost equally strong peaks at 3.2 meV and 4.4 meV 
which points, in agreement with earlier neutron data, towards a Tg quartet ground state. Further 
possibilities like a quasi-quartet ground state are discussed. 

PACS numbers: 71.27.+a, 75.10.Dg, 78.70.Nx, 75.30.Mb 



I. INTRODUCTION 

YbInNi4 came into the focus of interest when the first 
order valence transition from a trivalent high tempera- 
ture to an intermediate low temperature phase was dis- 
covered in YbInCu4.— In many rare earth compounds the 
4f electrons couple with the conduction electrons, lead- 
ing to the wealth of properties from magnetic order, to 
Kondo and Heavy Fermion states with non-conventional 
superconductivity and/or non-Fermi liquid behaviour, 
and intermediate valency. Hence the knowledge of the 
low energy excitations, especially the ground state wave 
function of the 4f electrons involved, is crucial for any 
further understanding of what mechanism leads to what 
property. The crystal-field splitting of YbInCu4 is diffi- 
cult to study directly because of its intermediate valent 
state at low temperatures. Instead YbInNi4 has been 
investigated, because it is trivalent over the entire tem- 
perature range. However, there are many discrepancies 
between the different crystal-field proposals^— so that we 
re-investigated the crystal-field scheme of YbInNi4 with 
polarized neutron diffraction and inelastic neutron scat- 
tering. 

At high temperatures YbInNi4 shows Curie- Weiss be- 
haviour with the full Yb 3+ magnetic moment and orders 
magnetically at 3 YbInNi4 forms in the cubic C15b 
Laves structure with T^ symmetry at the Yb site, so 
that the 8-fold degenerate Hund's rule ground state of 
Yb 3+ with 3=7/2 and J a =| ± 7/2), | ± 5/2), | ± 3/2), 
| ± 1/2) is lifted by the crystal- field into two Kramer's 
doublets \Tq) and ^7), and one quartet \Tg). In cu- 
bic site symmetry the J z admixture of these crystal-field 



wave functions is fixed so that the crystal-field scheme is 
fully determined with the crystal-field transition energies 
and the sequence of states. In Stevens approximation the 
transition energies and sequence of states is given by the 
Stevens parameters. For Yb 3+ with cubic T^ symmetry 
the two Stevens parameters and B§ fully describe the 
crystal-field scheme! 7 - 

There has been a controversy about the crystal-field 
schemes for several years and it became of interest again 
because of the more recent discovery that applying pres- 
sure (2.45 GPa) to YbInCu4 suppresses the valence tran- 
sition and leads to magnetic order at 2.4 Ki^r— There are 
several crystal-field propositions around which are briefly 
summarized in the following: 1) Inelastic neutron scatter- 
ing data by Severing et al£ showed a double peak struc- 
ture at low temperatures and the excitations at about 3 
and 4 meV were interpreted as transitions from a quar- 
tet ground state because the transition matrix element 
between the two doublets is zero^ 2) Sarrao et al. con- 
cluded from their entropy findings in the specific heat 
that the ground state should be a doublet and they spec- 
ified from magnetization data that it should be the ^7) 
state. Best fits were obtained when assuming the second 
doublet at 5.5 meV, and the |r 8 ) quartet at 10.5 meV£ 
3) From rare earth doped LuInNi4 data of the magnetic 
susceptibility and from ESR measurements Pagliuso et 
al. suggest also a ^7) ground state with the |Tg) quartet 
at 4 meV and the \Tq) doublet at 9 meV— . 4) Aviani et al. 
and Park et al. investigated the crystal-field ground state 
of YbInCu4 , where Aviani et al. suppressed the valence 
transition with 50% Y doping and Park et al. by ap- 
plying pressure. Aviani et al describe their specific heat 
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data well with a quasi-quartet ground state which is the 
inversed scheme of Severing et at, while Park et al. de- 
scribe their specific heat data under pressure best with 
the \Tg) quartet scenario and a total splitting of 28 K 
when taking the Kondo effect into account£^ which is 
in agreement with the neutron proposal for YbInNi.4. 



II. EXPERIMENTAL AND ANALYSIS 

These different proposition challenged us to make an- 
other attempt to determine the crystal-field scheme of 
YbInNi4 with neutron scattering. We have measured the 
low temperature magnetic form factor, which probes di- 
rectly the ground state wave function as Fourier trans- 
form of the spatial distribution of the 4f electron mo- 
ment, i.e. it is directly sensitive to the anisotropy of the 
crystal-field ground state. Moreover, we have performed 
inelastic neutrons scattering experiments on a powder 
sample with emphasis on measuring excitations at low 
temperatures with a resolution better than in reference^. 
Crystals were grown by flux growth and the CI 5b struc- 
ture has been verified by powder x-ray diffraction. 



A. Magnetic Form Factor 

The magnetic scattering intensity in an elastic neu- 
tron scattering experiment is determined by the magnetic 
structure factor Fm where 



F m = ^Y1 \m n \(Q)exp(iQR n )exp(-W„) (1) 



and where m(Q) is the Fourier transform of the mag- 
netization density at the scattering vector Q. The sum- 
mation is over the n atoms of the magnetic unit cell, 
exp(—W n ) is the Debye-Waller factor and R n the po- 
sition of the nth atom, [Ib is the Bohr magneton, and 
2 

r o = 7^5- with 7=-1.92 the gyromagnetic ratio of the 
neutron. In dipole approximation m(Q) depends only 
on |Q| and can be written as the product of the mag- 
netic moment m and the spherical magnetic form fac- 
tor /(IQI), m(|Q|) = m/(|Q|). /(|Q|) is listed in text 
bookSi 1 — However, the dipole approximation is not valid 
for large momentum transfers and/or when the spatial 
distribution of magnetic moments is strongly anisotropic 
as for 4/ moments in the presence of a crystal-field. 
Then the magnetic scattering intensity is proportional 
to |mj_(Q)| 2 , which is the square of the projection of 
m(Q) perpendicular to Q. This introduces in addition 
to the modulus a vector Q dependence to the magnetic 
intensities which does not exist in dipole approximation. 
This vector Q dependence of mj_(Q) can be used to de- 
termine crystal-field ground state wave functions (see e.g. 
reference UJ-UJ). Here we determined the vector Q de- 
pendence of m(Q) of Yb 3+ with the aim to obtain the 
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FIG. 1: (color online) Ytterbium magnetic amplitudes jm(Q) 
measured with an applied field of 2 T at T = 5 K (a) and T 
= 20 K (b). The full blue squares correspond to (001), the 
full red circles to (hhO), and the black open squares to other 
reflections with even indices. 



ground state wave function of YbInNi4. A polarized neu- 
tron diffraction experiment was performed on a YbInNi4 
single crystal using the D3 diffractometer at the high flux 
reactor of the Institut Laue-Langevin in Grenoble. The 
peak intensity of the Bragg reflections which is measured 
for neutrons polarized parallel and antiparallel to the ap- 
plied magnetic field, leads to the so-called flipping ratio 
R = I + 1 I~ . To deduce the magnetic contribution from 
R, a good knowledge of the nuclear structure factors i<N 
is necessary. This was achieved by an experiment on the 
same crystal on the 4-circle diffractometer D9 at ILL. 

The nuclear intensity measurements with unpolarized 
neutrons (D9) were performed at T = 5 K with two differ- 
ent wavelengths A = 0.84 A(II46 Bragg reflections) and 
A = 0.51 A(698 Bragg reflections). The crystal structure 
of our sample was refined to the C15b structure with 
the program MXDJI The scattering lengths were taken 
from the BNL tables^ as 12.43 fm for ytterbium and 
10.3 fm for nickel. As the scattering length of indium 
depends strongly on the wavelength, its real and imag- 
inary parts (b and respectively) were calculated for 
each wavelength from the values given at 1.8 A. The to- 
tal absorption coefficient \i was also calculated and all 
these wavelength dependent values are gathered in ta- 
ble I. The nickel atomic position XNi, the Debye-Waller 
thermal parameters W n , and the extinction parameters, 
block size t and mosaicity were refined, leading to 
consistent values at both wavelengths. These values are 
given in table II. 

The polarized neutron diffraction experiments (D3) 
were performed at two temperatures, 5 and 20 K, in a 
field of H = 2 T applied along the [1 10] axis of the crys- 
tal. The selected wavelength was A =0.825 A, and two 
erbium filters were used to suppress the A/2 contamina- 
tion. To optimize the extinction corrections ) 20 ' 21 flipping 
ratios of some particular reflections were also measured 
at two lower wavelengths A = 0.74 A and A = 0.52 A, 
thus taking advantage of the hot source available on D3. 
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TABLE I: Values of the indium scattering length (real part 
bin and imaginary part 6"i n ) and of the total absorption co- 
efficient £iYbinNi 4 as a function of the wavelength A. 



Device 


A 


bin 


b"ln 


MYbInNi 4 




(A) 


(fm) 


(fm) 


(cm- 1 ) 


D9 


0.84 


0.3924 


-0.0060 


1.42 




0.51 


0.3537 


-0.0083 


1.21 


D3 


0.825 


0.3917 


-0.0061 


1.42 




0.74 


0.3867 


-0.0063 


1.34 




0.52 


0.3564 


-0.0081 


1.21 



In YbInNi4, the nuclear and magnetic structure fac- 
tors are complex. Their imaginary part has two origins: 
the imaginary scattering length of indium and the non- 
centrosymmetric structure of YbInNi4. This second fea- 
ture affects the ytterbium contribution to the structure 
factor only for Bragg reflections with h,k,l odd since the 
Yb occupy the 4c lattice sites at Note that the 

expectation value mj_(Q) is real for any 4f wave func- 
tion, so that the magnetic structure factor Fm is real for 
h, k, I even. For m(Q) || H the nipping ratio R can then 
be written (see e.g. Ref. I22T ): 

R = iff + ^n 2 + 2 sin 2 aF^Fu + sin 2 aF^ 
F S + F n 2 ~ 2 sin2 vFnFm + sin 2 aF^ 

Fjtj and F^ are the real and imaginary parts of the 
nuclear structure factor F^ and a is the angle between 
the magnetic moment m(Q) and the scattering vector 
Q. This formula undergoes appropriate corrections in or- 
der to take into account instrumental imperfections (po- 
larization of the incident beam, nipping efficiency) and 
extinction effects. 

55 non-equivalent reflections were measured at T = 
5 K and 49 ones at T = 20 K, up to sin6»/A = 0.94 A -1 . 
Among them, 13 independent reflections measured at 
several different wavelengths were used to refine the ex- 
tinction parameters. With fixing the Ni positions Xjsn, 
the Debye Waller parameters W n , and block size t to the 
values deduced from the refinements on D9, the mosaic- 
ity g was refined to a value in complete agreement with 
the previous ones (see table II). The final values of the 
different parameters used for the data analysis are also 
given in that table II (see reference [20lf2ll ) . 

The Fm values deduced from the measured flipping 
ratios are directly related to the Fourier transform of the 
magnetization density m(Q) for the scattering vector Q 
(see eq.(l)). Fig. 1 shows the obtained |m(Q)| values 
at 5 K and 20 K versus sin0/X with |Q|=47rsin#/A. At 
5 K there is an anisotropy, i.e. all Bragg reflections lie 
between the (hhO) and the (OOZ)-type ones. The (00Z) 
reflection decrease steeper with sin 0/ A than the (hhO) 
ones. At 20 K, this anisotropy has nearly vanished. 

Only these special reflections are drawn in Fig. 2. 
We have now simulated |m(Q)| for the different crystal- 
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FIG. 2: (color online) Comparison of simulated ytterbium 
magnetic amplitudes |m(Q)| at T = 5 K for different ground 
states. The full symbols are the measured values for (00Z) 
(blue squares) and (hhO) (red circles), the lines are simula- 
tions for the (00Z) (blue) and (hhO) (red) directions for the 
different crystal-field propositions. The black line corresponds 
to the spherical form factor. The values for sin(9/A = (black 
triangles) come from magnetization. 

field scenarios as suggested by£~— The simulation was 
performed with the program package McPhase 2 ^ which 
includes form factor calculations beyond dipole approxi- 
mation. Fig. 2 displays the measured |m(Q)| values for 
the (00Z) and (hhO) reflections as full symbols. The value 
for sm0/X = is taken from magnetization. The solid 
lines are simulations. For the current plot the simula- 
tions have been scaled to the value at sin#/A = 0. The 
scaling factors are 1.07, 1.05, and 0.95 for the |r 6 ), |r 7 ), 
and |Tg) simulation. The simulations show that a ^7) 
ground state would yield the wrong anisotropy: For a 
\Tf) the (hhO) reflections decrease steeper with sm0/X 
than the (00Z) which is in contradiction to our observa- 
tion. Hence our form factor measurements clearly rule 
out the possibility of \Tr) as ground state. In contrast, 
the simulations for the \Tq) ground state of the quasi- 
quartet scenario and the for the \Tg) quartet both give 
the correct anisotropy. 



B. Inelastic Neutron Scattering 

Inelastic neutron scattering experiments on polycrys- 
talline samples are the most common technique to deter- 
mine crystal-field excitations in rare earth compounds. 
Here we present inelastic neutron data of polycrystallinc 
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TABLE II: Results of the refinements of the nuclear structure. Values of the nickel atomic position x^i, the Debye-Waller 
thermal parameters W n , extinction parameters t and g, and ~x 2 = Pi(A° bs — A\ a ) 2 /N t s — N var ) with Pi — 1/cr 2 and A t = I t 
(intensity) for D9 or Ai — Ri (flipping ratio) for D3. 
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D 9- refined 


0.84 
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FIG. 3: (color online) Scattering function of YbInNi4 for Q 
in between 0.35 A -1 and 2.3 A -1 at 2 K (open black cir- 
cles) and 5 K (full blue circles), and of LuAuCu4 (open red 
squares), also at 5 K. The insets represent the elastic scat- 
tering divided by 42. The extra scattering in the YbInN4 
data between 2.5 meV and 6 meV is attributed to magnetic 
scattering. 



YbInNi4 which were taken at the time-of-flight spectrom- 
eter MARI at the pulsed neutron source ISIS with an 
incident energy of 12 meV and a resolution of 0.4 meV 
(FWHM) at elastic scattering. Detectors from 29 = 8° 
to 52° are grouped together, resulting in an averaged Q 
vector at elastic scattering of 1.34 A. The sample was 
mounted in a cryostat and data were taken in the mag- 
netically ordered phase at 2 K and in the paramagnetic 
phase at 5 K, 20 K, and 40 K. 

Fig. 3 shows the scattering function S(Q,uj) versus 
energy transfer for YbInNi4 at 2 and 5 K, and of the 
non-magnetic isostructural compound LuAuCu4, also at 
5 K. The data are corrected for detector efficiency, ab- 
sorption, and are scaled to sample mass so that the scat- 
tering intensities are comparable. Due to the smaller nu- 
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FIG. 4: (color online) Temperature evolution of the scattering 
function of YbInNi4 for 12 meV incident energy and |Q| in be- 
tween 0.35 A -1 and 2.3 A -1 . The insets represent the elastic 
scattering divided by 42, the dashed lines crystal-field excita- 
tions resulting from a crystal-field fit, and the solid lines are 
the total fitted magnetic scattering. Note, the fitting region 
was restricted to energy transfers below -2 meV and above 
2.5 meV to avoid the strong elastic and the spurious scatter- 
ing at the tail of the elastic line. 
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clear cross-section of LuAuCu4 with respect to YbInNi4 
the elastic scattering of the latter is stronger. The Lu 
data prove that phonon scattering is negligible in the 
present energy window and detector grouping. The ex- 
tra scattering of the YbInNi4 sample between 2.5 meV 
and 6 meV is therefore identified as magnetic scattering. 
The shoulder which appears on the neutron energy loss 
side of the elastic line appears in the spectra of the mag- 
netic and non-magnetic sample, and is attributed to spu- 
rious Bragg scattering in either sample which has been 
reflected from the cryostat walls, reaching the detectors 
time delayed, and therefore in an inelastic channel. Since 
it is difficult to account for this quantitatively and since 
we know from previous high resolution data£ that the 
quasielastic scattering of YbInNi4 is too narrow to be re- 
solved in the present experiment we concentrate on the 
inelastic scattering above 2.5 meV and on the neutron 
energy gain side for T = 20 K and 40 K. 

The magnetic scattering intensity between 2.5 meV 
and 6 meV consists clearly of two lines. At 5 K these 
lines can be fitted with two Lorentzians centered at about 
3.2±0.f meV and 4.4±0.1 meV. In the magnetically or- 
dered phase at 2 K the lower crystal-field excitation ap- 
pears at slightly larger energy transfers with respect to 
5 K which is probably due to the influence of the mag- 
netic order at 2 K. The line widths remain unchanged. 

Fig. 4 shows the temperature evolution of the scatter- 
ing function in the paramagnetic phase from 5 K, 20 K 
to 40 K. The double peak structure survives up to 40 K 
without a shift in energy, although the excitations be- 
come broader as temperature rises. The present incoming 
neutron energy of 12 meV provides an energy window up 
to 10 meV. At none of the temperatures magnetic scat- 
tering was detected at energy transfers larger than than 
4.4 meV. 

The scattering function S(Q,w) can exhibit two mag- 
netic excitations at 2 and 5 K only if 

a) the ground state is the \Tg) quartet and two ground 
state excitations into the \T§) and \Tj) doublets take 
place, or 

b) the ground state is a quasi-quartet and the two dou- 
blets which are 1 meV apart are sufficiently populated, 
or 

c) if the crystal-field states are Zeeman split due to the 
proximity of the magnetically ordered state. Then even 
a doublet could give rise to two ground state excitations. 

d) Or if some structural distortion is present, so that the 
site symmetry is no longer cubic. 

Proposition a) is based on the fact that the transi- 
tion matrix element between the two doublets is zero. 11 
This new set of data is also, like the previously reported 
neutron data, well described with a \T$) quartet ground 
state. For a quantitative analysis the magnetic form fac- 
tor has to be taken into account. Here, for the purpose of 
describing non dispersive excitations in a powder sample 
it is sufficient use the spherical, text book magnetic form 
factor /(IQI). 1 - The result of a quantitative crystal- field 
analysis where all three temperatures were fitted simulta- 
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FIG. 5: (color online) Simulated inelastic scattering function 
at 5 K for a quasi-quartet (left) and quartet (right) ground 
state crystal-field model. 

neously is shown in Fig. 4 (see black lines). For both exci- 
tations the same line width was assumed. The fitted line 
widths increase with temperature from 0.54±0.05 meV 
at 5 K to 0.7±0.05 meV at 20 K and 0.9±0.08 meV at 
40 K. The best fit was obtained for B° = -1.110 -3 meV 
and Bg = — 8.8810~ 5 meV, which is consistent with ref- 
erence |3| but seemingly in conflict with the entropy find- 
ings of the specific heat 2 . Although b), the quasi-quartet 
scenario a la Aviani et ai£ is tempting, but a simple 
simulation of the scattering function for a quasi-quartet 
ground state excludes this possibility as an explanation 
for the inelastic neutron data of YbInNi4 : neither at 2 
nor at 5 K the first excited doublet at 1 meV is sufficiently 
populated to give rise to two almost equally strong exci- 
tations. This is demonstrated in Fig. 5 where the result 
of a crystal-field simulation with a quasi-quartet ground 
state is shown for 5 K. Possibility c) of a molecular field 
split doublet ground state could give rise to two excita- 
tions at low temperatures, but for e.g. a \Tq) ground 
state the same molecular field of «0.33 meV would have 
to be assumed for the magnetically ordered and the para- 
magnetic phase. Another difficulty with this assumption 
is that at higher temperatures, when no Zeeman split- 
ting and intermixing of states due to a molecular field 
is present, the spectra should be shifted in energy with 
respect to the low temperature data. This has not been 
observed. Possibility d) of a structural distortion can be 
ruled out from our structural analysis which is part of 
the form factor measurement and from X-ray data which 
were taken in order to verify sample quality. Hence, this 
new set of neutron data confirms the previous inleastic 
neutron scattering results of a |r 8 ) quartet ground state. 



III. DISCUSSION 

The observed anisotropy of the magnetic form factor 
is clearly not compatible with the spatial distribution 
of a |r 7 ) ground state although magnetization and sus- 
ceptibility measurements^ favour ^7). This demon- 
strates the difficulty to determine crystal-field ground 
states from magnetic measurements where thermal aver- 
aging, exchange, possible Kondo interactions, etc. enter 
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the description but are not easily taken care of. The mag- 
netic form factor, however, cannot distinguish between 
r 6 ) and |r 8 ) because the projection of the orbitals of 
these states on to the scattering plane is very much the 
same. A \Tq) ground state would be in agreement with 
the entropy findings of the specific heali^ but as has been 
pointed out in the upper paragraph, it fails to explain the 
double peak structure in the inelastic neutron scattering 
data. All inelastic neutron scattering data exhibit a dou- 
ble peak structure at all temperatures and can be ana- 
lyzed consistently with a \Tg) quartet ground state and 
two ground state excitations. From the present experi- 
ment with 12 meV incident energy and from the previous 
one with 17 meV we can conclude further that there is 
no additional ground state excitation within an energy 
window up to 15 meV. 

However, it should be mentioned that additional peaks 
in the inelastic neutron data have been observed for ex- 
ample in CeAl2 ~ Here strong electron phonon coupling 
gives rise to an additional peak in S(Q,w). But this usu- 
ally happens at higher energy transfers of acoustic zone 
boundary or optical phonon branches because the obser- 
vation of such a bound state requires a certain phonon 
density of states which is not given in the range where 
acoustic branches rise steeply. For CeAl 2 the bound state 
is observed at 10 meV which coincides with the first peak 
in the phonon density of states. Ytterbium is not much 
heavier than cerium so that it is unlikely to have a fiat 
phonon branch in the energy range of 3-4 meV. Hence we 
consider it unlikely that electron phonon coupling gives 
rise to an extra excitation. 

We consider the quartet ground state in YbInNi4 
as the most likely solution. This makes us speculate 
whether YbInNi4 is a candidate for quadrupolar order 
as e.g. in CeBg. We further believe that the mag- 



netic exchange interaction in YbInNi4 is predominantly 
antiferromagnetio^ and not ferromagnetic as suggested 
in reference which then leads us to speculate whether 
there is magnetic frustration due to the fee crystal- 
structure of YbInNi4 where the Yb ions are arranged 
in corner sharing tetrahedra. We therefore do not think 
that the entropy findings of Rln2 by Sarrao et al. arc 
necessarily in conflict with a quartet ground state. 



IV. SUMMARY 

The crystal-field scheme of YbInNi4 has been deter- 
mined with magnetic form factor and inelastic neutron 
scattering measurements. The asymmetry of the mag- 
netic form factor clearly rules out a ^7} doublet as 
ground state (see Fig. 2) and would be in line with ei- 
ther a |r@) or \Tg) ground state. A quasi-quartet ground 
state, where the two doublets are 1 meV apart, would not 
yield two almost equally strong excitations in the inelas- 
tic neutron data at 5 K since the first excited state would 
not be sufficiently populated (see Fig. 5). The inelastic 
neutron data confirm the existence of two crystal-field ex- 
citations at 3.2 and 4.4 meV at low temperatures and due 
to selection rules this is only possible for a |Tg) quartet 
ground state. 
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